DTP/96/22 
RAL-TR-96-019 
March 1996 



yQ-, Can partons describe the CDF jet data? 



oo 



> 

m 
m 
o 

a,: 

Oh! 



E. W. N. Glover", A. D. Martin", R. G. Roberts 6 and W. J. Stirling a ' c 

a Department of Physics, University of Durham, Durham, DH1 3LE 
b Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX 
Department of Mathematical Sciences, University of Durham, Durham, DH1 3LE 



Abstract 
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^ ■ The recent CDF single jet inclusive measurements at Fermilab are incorporated in a 

global next-to-leading order parton analysis of the available deep inelastic and related 
data. We find that it is impossible to achieve a simultaneous QCD description of both 
the CDF jet distribution for transverse energies Ej> > 200 GeV and the deep inelastic 
structure function data for x > 0.3. However, the CDF data for Et < 200 GeV and the 
deep inelastic data are adequately described provided that the QCD coupling as(Mg) is 
increased from its preferred deep inelastic value to as(M§) ~ 0.116 — 0.120. 



The measurement of the differential cross section for inclusive central jet production at the 
Fermilab pp collider has recently been reported by the CDF collaboration M for jet transverse 
energies, Et, in the range 15 to 440 GeV. At the higher Et values these measurements probe 
the substructure of the proton in a previously unexplored kinematic region, equivalent to 4- 
momentum transfer squared Q 2 ~ 10 5 GeV 2 which corresponds to distance scales of O(10~ 17 m) 
or less. Intriguingly, the experimental results for Et > 200 GeV show evidence of a possible 
deviation above the behaviour predicted by next-to-leading order (NLO) QCD based on the 
current sets of parton distributions — distributions which are obtained from global analyses 
|fj, ||, f| of a wide range of deep inelastic and related data. Clearly before explanations based 
on New Physics can be taken seriously, it is crucial to see if the parton distributions can be 
adjusted to accommodate the jet measurements whilst retaining a satisfactory description of 
the other data. Here we address this vital question. 

We begin by comparing the CDF jet measurements with the (NLO) prediction obtained 
from the MRS (A') set of partons |3J which, at present, seem best able to describe all the 
other data. The calculation uses the next-to-leading-order parton level Monte Carlo JETRAD 
and the cuts and jet algorithm applied directly to the partons are modelled as closely as 
possible to the experimental set-up. The jets are defined according to the Snowmass algorithm 
with a jet cone size R = 0.7 0, and are required to lie in the pseudo-rapidity range between 
0.1 < \rj\ < 0.7. The factorization and renormalization scales are chosen to be \ip = = Et/2. 
The fractional difference between data and theory is given by comparing the data points with 
the horizontal line in Fig. 1.Q Over the observed range of Et the experimental cross section 
falls by more than 8 orders of magnitude and the sum of the correlated systematic uncertainties 
grows from ± 18% at E T = 50 GeV to ± 28% at E T = 300 GeV. However, the statistical 
precision of the data is such that the difference between the observed and predicted shapes 
of the Et distribution must be studied seriously. Although we have selected MRS (A'), it is 
important to note that the comparison is very similar for the other parton sets obtained from 
global analyses which include both fixed-target and HERA deep inelastic data, see Fig. 1 of 
Ref. |!]]. Disregarding a small change in normalization, the nature of the discrepancy between 
the MRS(A')-based prediction and the data differs according to whether the jet E T is below 
or above about 200 GeV. Interestingly it is around this value that jet production changes from 
being dominated by gg-initiated to (valence quark) gg-initiated QCD subprocesses, see Fig. 2. 

Before we present a global parton analysis which incorporates both the deep inelastic and 
jet data we discuss some general trends that can already be established, bearing in mind that 
a centrally produced jet of transverse energy Et samples partons at 



where -^i =1.8 TeV at the Tevatron collider. First we note that the discrepancy in the shape 
of the jet distribution for Et < 200 GeV could be removed either by significantly increasing the 

^or clarity we have not made the small overall renormalization of the A' prediction which would improve 
the average absolute value, but not the description of the shape of the jet distribution. 
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QCD coupling from the MRS (A') value of as(M§) = 0.113 (see Refs. || or by increasing 
the gluon in the x ~ 0.1 region (see, for example, the MRS(Dq) [nj "base" line in Fig. 1 of 
Rcf. ||T|). Both options lead to disagreement with data elsewhere. The effect of increasing as 
is to increase the rate of increase in the partons at small x and to steepen the low Ex jet cross 
section as required, but it also leads to scaling violations of deep inelastic data which are too 
rapid. The second option is contrary to HERA data which require the gluon to be larger (and 
steeper) than D' for x < 10 -2 and hence, by momentum conservation, smaller (and steeper) 
for x ~ 0.1, see Fig. 15 of Ref. @. 

Turning now to the discrepancy for > 200 GeV, we note that in this region the jet 
distribution depends on the partons (mainly the valence quarks, but also the gluon) with 
x ~ 0.35 — 0.5 and Q 2 ~ 10 4 — 10 5 GeV 2 . The discrepancy has been accommodated in a recent 
analysis Jl I J by modifying the gluon distribution in this x region. Traditionally the gluon is 
constrained in the region 0.35 < x < 0.55 by the WA70 prompt photon data pf . In Ref. [|ll[ 



it is argued that the modified gluon can be made consistent with the WA70 data by taking 
into account scale uncertainties and kx broadening. However, increasing the gluon in this x 
range implies a decrease in the region x ~ 0.01 — 0.3 (see Fig. 2 of |TTJ) so the resulting partons 
do not describe the shape of the more precise jet data with E T < 200 GeV, see the dotted 
curve in Fig. 1. Alternatively, it might appear that a judicious increase of the valence quark 
distributions by about 10% in the region x ~ 2E T /y/s ~ 0.35 - 0.45 and Q 2 ~ 10 4 - 10 5 GeV 2 
would remove the discrepancy for Ex > 200 GeV. However, these distributions are constrained, 
through GLAP evolution, by deep inelastic measurements at lower Q 2 . In particular, precise 
large-a; structure function data from BCDMS (Ff,F^ d ) |TJ, CCFR (Fg N ,xF£ N ) 0J and 
NMC {F^ n / F^) Jl5| pin down the dominant u and d distributions very accurately, see for 
example Ref. @, with combined statistical and systematic errors of the order of a few percent. 
Data at smaller x also, through the effects of evolution and sum rules, constrain the structure 
functions at larger x. The situation is further complicated since an increase of the value as 
(which increases the jet subprocess cross section) tends, through evolution, to decrease the 
valence, and hence the jet, distributions at large Q 2 . 

It is clear that the jet data must be included in a global analysis to investigate whether 
a compromise overall fit can be achieved by adjusting the partons and as(Mf ), such that we 
have a satisfactory QCD description of the detailed shape of the entire jet spectrum without 
destroying the fit to the other data. We therefore perform a next-to-leading (NLO) analysis of 
the available deep-inelastic and related data as described in Ref. 0, but updated to include the 
published measurements of F2 at HERA [[L6| together with the CDF jet data [[[]. The inclusion 
of the HERA data impose an important constraint on the gluon at small x (and rule out, for 
example, partons such as MRS(Dq) |10| which were found to give an excellent prediction of the 
jet shape for Et < 200 GeV). The HERA data pin down the quarks and gluons at small x and 
thus (through the sum rules) restrict their variation at larger x. 

We fit to the shape of the jet data in the region Et > 50 GeV where the NLO QCD descrip- 
tion is, to a good approximation, scale independent (and where the experimental ambiguities 
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from the underlying event are small [II]). To implement an economical NLO description of the 
inclusive jet distribution within the global fit, we first calculate the K factor, 

Tyr &LO + &NLO 
K = 

0"LO 

as a function of E? using MRS(A') partons, but with larger as(M§) = 0.12. The renormal- 
ization and factorization scales are set equal to Et/2. A global analysis is performed in which 
the jet data are fitted to the distribution Ka^o- The resulting optimum set of partons to- 
gether with the corresponding value of as(M§) are then used to recompute the K factor, and 
the global fit repeated. This iterative procedure converges rapidly. For the above choice of 
scales the K factor is found, to a good approximation, to be independent of E T in the interval 
50 < Et < 400 GeV and equal to 1.10. The whole analysis is repeated taking the renormal- 
ization and factorization scales equal to Et- It is found to yield essentially the same global fit, 
but the normalization of the jet data increases. 

The new set of partonsf| which we denote by J, is found to have a significantly larger 
value of the QCD coupling, as(M§) = 0.120, as compared to 0.113 of MRS(A'), in order to 
give an improved description of the precise CDF jet data with E T < 200 GeV, see Fig. 1. 
The improvement is obtained at the expense of a somewhat less satisfactory description of the 
scaling violations of the fixed-target deep-inelastic data. This is evident from the \ 2 values 
listed in Table 1. As an example, Fig. 3 shows the description of the large x BCDMS deep 
inelastic data JT3|. (A similar plot exists for the F^ and F% N ,xF^ N data.) Note that there is 
in fact a continuum of such fits, spanned by A' and J, in which a s (M§) increases from 0.113 to 
0.120 for which the fits to the deep-inelastic (low E T jet data) gradually deteriorate (improve). 
In selecting fit J we are deliberately giving more weight in the fit to the CDF data. 

The compromise fit, J, does not improve the description of the CDF jet data with Et > 
200 GeV. In fact we see from Fig. 1 that the fit is marginally inferior to A'. It turns out to be 
impossible to obtain a reasonable description simultaneously of the deep inelastic and CDF jet 
data over the full range of Et- Because of their smaller errors the low Et CDF data points have 
a bigger influence on the fit than those at large Et- To illustrate the conflicting requirements of 
these data we successively omit deep inelastic dataQ starting from those at the largest x values 
until we can achieve a satisfactory description of the entire CDF jet Et spectrum. We find 
that it is necessary to omit data with x > 0.2, that is a major fraction of the deep inelastic 
data, before a reasonable fit to the jet spectrum is obtained. The partons obtained in this way 
are denoted by J'. Although the fit to the jet Et spectrum is satisfactory (see Fig. 1), these 
partons give a poor description of the omitted deep inelastic data — see, for example, the J' 
prediction for the BCDMS data shown in Fig. 3. 

In the J' fit the quarks are drastically modified at large x compared to MRS (A') in order 
to accommodate the large Et CDF jet data. The question addressed in Ref. |IlJ is whether 

2 A similar fit, but with the jet data described by ulo was obtained in Ref. |L7| . 
3 However, we keep the NMC measurements [|l5| of FJ 1 /Ff in the fit. 
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Table 1: % 2 values for various subsets of data; the values shown in square brackets correspond 
to data that are omitted from the fit. 
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the quarks can be kept essentially unchanged and the large E T excess explained solely by a 
modified gluon density at large x. We have also performed fits of this type and find that the E T 
distribution cannot be adequately described below and above E T ~ 200 GeV simultaneously - 
in qualitative agreement with the conclusion one draws from the dotted curve [TT] of Fig. 1. We 
have also explored the possibiliity of modifying the scale for as by including a factor (1 —x)/x, 
which helps give the jet cross section a rise with increasing Et- However, the resulting fit can 
produce qualitative agreement only over a very limited region (E T > 200 GeV) provided a large 
renormalization of the jet data is made. 

Fig. 4 compares the quark and gluon distributions, at Q 2 = 10, 10 4 GeV 2 , of the new sets 
with those of MRS (A'). We plot the singlet quark distribution, J2i(<li + Qi), since this is most 
closely related to the combination which is sampled by the jet cross section. Note that the area 
under the curves is the total momentum fraction carried by the quarks and gluons. We see 
that the main difference between the A' and J distributions comes from the different value of 
as- the evolution between Q 2 = 10 GeV 2 and 10 4 GeV 2 is more pronounced for the J partons, 
such that at the latter scale the J quarks are respectively smaller and larger than the A' quarks 
for x > 0.1 and x < 0.1. This difference leads to a better fit to the small Ej- CDF data, see 
Fig. 1. The starting gluons are very similar in the two sets. For both the A' and J partons 
the overall renormalization required to give the best fit to the CDF jet data is well inside the 
experimental uncertainty. The J' partons, on the other hand, have significantly harder quark 
distributions at large x, as anticipated from Figs. 1 and 3. At Q 2 = 10 4 GeV 2 , the J' quarks 
are larger than the A' quarks for x j> 0.35 and x < 0.08, giving an improved fit to the CDF jet 
data both at large Et and at small Et- Notice also that the J' gluon is significantly smaller 
than the A' gluon at small x and Q 2 = 10 GeV 2 . This is required in order to maintain a good 
fit to the HERA measurements of dF 2 /dlnQ 2 ~ asg with a significantly larger as- 

In summary, we have attempted to incorporate the CDF inclusive jet cross section data in 
a global parton distribution analysis. We find that it is impossible to accommodate both the 
jet data over the complete Et range and the deep inelastic structure function data.0 However, 
the CDF data for Et < 200 GeV and the deep inelastic data can be reasonably well fitted 
simultaneously. The main effect of including the jet data is to increase the value of as from 
the standard deep inelastic value. Interestingly, this slightly improves the fit to the HERA and 
low-x fixed target data. The fit to the large-x fixed data is however considerably worsened, see 
Table 1. The result of the combined fit is the J set of partons. It is impossible to simultaneously 
fit the large Et > 200 GeV jet data. If one gives up the fit to the small Et data, one can 
achieve a reasonble large Et fit by adjusting the gluon distribution, as in the analysis of 
Ref. [TT| for example. A fit over the whole jet Et range (resulting in the J' set of partons) 
yields quarks which are completely incompatible with the large- x structure function data, see 



4 This conclusion assumes that the systematic errors of the CDF data are strongly correlated and that the 
statistical errors define the shape of the distribution. It is relevant to note that for the UA2 inclusive jet data jl8| 
there is no evidence for a discrepancy with the shape of the predicted NLO QCD distribution (i.e. a "break" in 
the region of xt = ^Et/^/s 0.25), although in this case the large combined statistical and systematic errors 
may mask such structure; see Fig. 5 of Ref. pjl. 
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Fig. 3, a value for as which is several standard deviations larger than the current world average, 
and a renormalization of the CDF jet data which is on the edge of the allowed range PJ. 

We therefore conclude that it is unlikely that the difference between the CDF inclusive jet 
cross section data and the standard NLO QCD prediction can be attributed to a deficiency in 
our knowledge of parton distributions. 
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Figure Captions 



Fig. 1 The fractional difference between the measured CDF inclusive jet cross section |]T| and 
the NLO QCD predictions based on MRS(A') M, J and J' parton sets. Also shown (by 



the dotted curve) is the gluon-modified description of Ref. ||11|| . The numbers shown in 
brackets are the renormalization factors used for the J and J' partons. Only the statistical 
errors of the data are shown. 

Fig. 2 The fraction of the leading-order jet cross section d 2 a / ' dETdr)\ v =Q at the Tevatron origi- 
nating from gg, qg and qq {+qq) initiated QCD subprocesses using the MRS (A') set of 
partons Ref. j4j. 

Fig. 3 The description of the large x BCDMS measurements of F% by the MRS (A') partons of 
Ref. pi, and the J and J' partons of this work. 

Fig. 4 The parton distributions g and 2i((Zi + <?i) at Q 2 = 10 and 10 4 GeV 2 . 



S 



1 .4 



.2 - 



1 - 



< 



o 

M 

d> 

o 
a 

CD 



0.8 



0.6 



13 0.4 

a 
o 



o 



0.2 - 



Comparison with CDF jet data 



MRS(A') (not renorm.) 
J(/1.036) 
J' (/1.205) 

Huston et al. 0.85 Jet-Fit 



♦+ "7 



-0.2 



-0.4 



50 100 150 200 250 300 550 400 450 

E T (GeV) 



Fig. 1 




Ep (GeV) 



Fig. 2 



0.3 



0.25 
0.2 
0.1 5 



Comparison with BCDMS data x=0.35 



7 ep 



0.2 
,0.18 



0.1 6 
0.14 
0.12 
0.1 



F? 



0.12 
0.1 
0.08 
0.06 




10 



Q 2 (GeV 2 ) 



10 





■v. 


x=0.45 




I 


III I 





* • • 

i ~ ~ - 



10 



Q 2 (GeV 2 ) 







x=0.55 






1 Q 2 (GeV 2 ) 


10 2 




0.06 
0.05 








= ^ x=0.65 " 






0.04 








0.03 








0.02 


_ i i iii 


ii 


1 ~ ' 


1 Q 2 (GeV 21 








10 2 





Fig. 3 




Fig. 4 



